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Abstract 

In this paper we are discussing the problem of low-energy greybody factors for fermions 
emitted by a Schwarzschild-de Sitter black hole. In our study we are using the analyti¬ 
cal methods proposed by Unruh some time ago for determining the greybody factors. We 
have found that at low energies the greybody factors are constant for a given total angular 
momentum (similar to what happens in the case of scalar particles reported before in the 
literature). Also our results are indicating an enhancement in the energy spectrum if one is 
increasing the value of the cosmological constant. These results are consistent with numer¬ 
ical calculations performed in Ref.[14] by S.F. Wu et al. (Phys.Rev.D 78, 084010). 
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1 Introduction 

Since the discovery more than a decade ago that our Universe has entered a new phase of 
accelerated expansion [1, 2] the study of non-flat asymptotically space-times geometries has 
drawn more and more attention from the scientific community. One of the most studied non¬ 
flat geometries is the de Sitter one due to it’s rich symmetries and the fact that if our Universe 
continues it’s accelerated expansion (which could be driven by the presence in the Einstein 
equations of a nonzero cosmological constant) then it will pass in the far future through a 
de Sitter phase. Moreover, the inflationary phase [3, 4] can be also described with a good 
approximation by a de Sitter geometry. It is worth noticing that Lagrangian modified gravity, 
which provides a dynamical mechanism for an effective cosmological constant implemented by 
the higher curvature terms in the Lagrangian, has found an application in cosmology (see e.g. 
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[5, 6]) in order to explain dark energy effects. Moreover, these so-called f{R) theories contain 
Schwarzschild-de Sitter metric as a vacuum solution [7]. They were also tested in the solar- 
system and galactic scales [8]. 

On the other hand the study of black holes also continues to be a hot topic in research today. 
Since Hawking’s [9] proposal that black holes can radiate energy through emission of particles, 
there have been a lot of research done in this line. However, it seems that the vast majority of 
studies concentrated on the asymptotically flat case while the non-flat case has been less studied 
so far. One of the quantities characteristic to the Hawking radiation are the greybody factors 
that describe the departure of the spectrum from the one of a pure black body. 

Thus in the present study we investigate the problem of low energy greybody factors for 
fermions emitted by a Schwarzschild-de Sitter (SdS) black hole. In order to do this we apply 
similar methods used by Unruh [10] for the case of Schwarzschild black holes (see also more 
recent studies [11, 12]). This consists mainly in obtaining approximative analytical solutions 
to the field equations in the vicinity of the black hole, respectively the cosmological horizons 
and then matching these solutions (technique valid only at low energies) and determining the 
absorbtion probabilities or equivalently the greybody factors. Given the complexity of the 
differential equations in the region between the two horizons we will perform our analysis only 
for small values of the cosmological constant. 

Up to our knowledge this is one of the first attempts to study analytically the emission 
of fermions from SdS black holes. Besides the numerical study performed in ref. [14] we 
haven’t found other studies that treat the problem of fermion emission in the 4-dimensional 
SdS geometry. The existing literature seems to be concentrated on the scalar case for which 
both numerical and analytical studies where performed [13, 14, 15, 19, 16, 17, 18, 20]; or on 
two-dimensional Vaidya-de Sitter toy models [21, 22, 23]; on thermodynamical radiation via 
tunneling [24] or the associated quasinormal frequencies [25, 26, 27]. 

The outline of the paper consists in presenting some general notions about the Dirac equation 
in curved space-times and a few basic proprieties of SdS black holes in section 2. The solutions 
to the field equations are presented in section 3 where we concentrate our attention especially 
on hnding the Dirac modes in an intermediate region between the two horizons present in SdS 
geometry. Section 4 is concentrated on finding the greybody factors for this geometry followed 
by a brief study of energy spectrum. We end our investigation with a brief section dedicated to 
our main conclusions. 

Trough this paper we use natural units with c = G = h = 1 and a (-|-,—,—,—) metric 
signature. 

2 Preliminaries 

In this section we present a short review of the Dirac equation written in the Cartesian gauge 
proposed some time ago in [28, 29]. The main advantage of this gauge is that it allows us to 
separate the spherical degrees of freedom in the same manner done for the central problems in a 
flat space-time [30]. This implies that the Dirac field will be a linear combination of particular 
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modes. Those of positive frequency and given energy e will read 


= a(r)[T+(r)$+ . ( 0 >) + F .(r)$ ( 0 >)]e 


—iet 


( 1 ) 


where we denoted by ^ the usual four-component angular spinors [30, 31] and by 
the radial part of the wave function; a{r) is a metric dependent function. The parameters mj 
and k in ( 1 ) stand for the standard quantum numbers and the eigenvalue k is related with the 
total angular momentum j and the orbital angular momentum I by the following forlmulas 




( 2 ) 


In particular, k = —1 corresponds to j = ^, / = 0 while k = 1 gives rise j = 1 = 1- 

The Dirac equation describing a free spinor field if: of mass m in a curved space-time back¬ 
ground can be defined in a local frame (characterised by the orthogonal tetrad-gauge fields 
and e“) as 

= 0, (3) 

The covariant derivative can be defined with the help of the spin connections 5"^ = |[ 7 “, 7 ^], 
where the 7 -matrices satisfy { 7 “, 7 ^} = 2 r/"^, 


04 = 40, = 4 + 

and where the local frame components of the connection read FT = a)> 

^ 1/3 t)eing the standard Christoffel symbols. 

The line element for a given manifold can be written in terms of the tetrad fields and 
the 1 -forms = eudx'^ as 

= 9)iudx^dx'' (5) 

Using the tetrads defined in ref. [28, 29], namely: 


a; 
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OJ 


OJ 
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dr 
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cos 9 

-;r^ sm 6 

u{r) 



v(r) 



rw{r) 

v{r) 

rw{r) 

v{r) 


sin 9 sin cj) dcj ), 
sin 9 cos (j) dcj ), 


( 6 ) 


the line element in a central static chart with spherical coordinates (t, r, 9, cp) can be expressed 
as 




w{r)‘^ 



dr'^ 

u{rY 



(dd^ -|- sin^ 9d(j)‘^) 


(7) 


where u(r), v{r) and w{r) are three arbitrary r-dependent functions. The form of a{r) intro¬ 
duced in ( 1 ) turns out to be a(r) = v{r)w{r)~^^‘^/r. 
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For a Schwarzschild-de Sitter space-time the line element is defined by 


ds^ = h{r) dt^ — —— — r‘^{d6‘^ + sin^ 9 dcfP '), 
h{r) 


( 8 ) 


with h{r) given by 


2 M A p 

hir) = 1- r 

^ ’ r 3 


(9) 


where M is the mass of the black hole and A stands for a positive cosmological constant. 
Comparing (7) with (8) we can make the following identifications 


u{r) = /i(r), v{r) = \/h{r) = w{r) 


( 10 ) 


It can be shown [28, 29] that the radial wave functions satisfy a radial eigenvalue problem 


HrJ- = eJ- written in a matrix form as 


m-/k(r) -/iWj; + f \/S(0 


+ -mv'TW 


riM 

KJ’-) 


r,+ (r) 


= e 


( 11 ) 




This system of first order differential equations is not analytically solvable due to the com¬ 
plicated form of h{r). However, in specific regions of the space-time we can find approximative 
analytical solutions. 

Before moving on to hnd solutions to eq. (11) let us say a few words about the properties 
of (8). Looking for the roots of h{r) which will point out the singularities of the metric (8) we 
find that it can have in general three different roots which will correspond to three horizons. 
However, only two of them are physical and they correspond to the black hole horizon (located 
at r = Tf,) and respectively to the cosmological horizon (located at r = Tc). The third solution to 
h{r) = 0 turns out to be related to the other two by rx = — (j'b + ^’c)- Given that the condition 
0 < ^ = AM^ < 1/9 is satisfied[32, 33] we will always have the two horizons satisfying r^, < rc- 
For the critical value of ^ = 1/9 the two horizons will merge and a Nariai[34] black hole will 
appear, while for ^ > 1/9 no horizons exist. In this paper we will study only the case of ^ < 1/9. 


3 Solutions to field equation 


The above coupled differential equations (11) can be transformed into a single second order 
equation. After making a change of variables defined by 


dr h 

dx 1 -|- A'v/h 


( 12 ) 


where A = m/e, the new equation for the upper component T“*'(r) reads 


d?F+ r 2 A - \Vh \ d ( kVh \ 
dx"^ ^ yi-FAvAy dx y(l \y/h)r j 


k^h 

(1 -|- \VhY 


F+ =0 


(13) 


In what follows we will solve this equation using the methods proposed by Unruh in [10] 


4 


















and developed further, for ex. in refs. [11, 12]. The solutions obviously depend on the value 
of the parameters: {k,M,A). However, in order to simplify the notation we shall skip writing 
them explicitly. 


3.1 Solutions near the two horizons and 

In the regions near the black hole horizon (located at r —)■ r;,), respectively the cosmological 

horizon (at r Vc) the function h{r) goes to zero so that Eq. (13) reduces to a more simple 

equation 

+ = 0 (14) 

which has the following general solution 

F+(x) = + (15) 


It is easy to see from (11) that the variable x behaves near the two horizons as x 
Integration constant can be chosen by imposing suitable initial conditions 


In h. 



2 A„ 


In/i = pln/i, if r 


n 


-1 


— lnh=qlnh, if r ^ Tc 


(16) 


Taking into account the ingoing boundary condition according to which only ingoing modes 
must exist in the vicinity of the black hole horizon, the solution in this (transition) region 
becomes 

F+ = « ^tr^-ieplnh 

At the cosmological horizon we have no restriction, thus the solution will be a combination 
of ingoing and outgoing modes 


C 


j^in ^—iex _|_ j^out^iex 

ieq In h _|_ j^out^ieqlnh 


(18) 


3.2 Solutions in the intermediate region ri, < r < Xc 


If we take into consideration the fact that in the intermediate region the terms proportional to 
or are much smaller than the other terms[10], then Eq. (13) for the radial wave function 
reduces to 


£f+ 

dx"^ 



^ kyfh ^ 

k'^h 

dx 

i^(l + \y/h)r J 

(1 + XyfhY 


F+ = 0 


(19) 


Our next goal will be to find an analytical solution to this equation. We start by introducing 
a new function H defined as 


dF~^ ky/h 
dx (1 + Xy/h)r 


( 20 ) 
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with the help of which Eq. (19) can be rewritten as 


k 


dH 


dr y/hr 


H = 0 


( 21 ) 


Due to the complicated form of h{r) Eq.(21) has no integrable solution. However, in the 
more simple case of Schwarzschild black hole for which 


h —y ho — 1 — 


2M 


( 22 ) 


equation (21) admits the following solution [10] 


Ha = C 


1 — y/hp 

1 + 


—k 


(23) 


Because the cosmological constant A has for the present expansion a very small value we 
can consider 6 = A/3 as a small parameter and expand Ij'/h around the value b=0 to obtain 


k 


k 1 /cr , 
+ --^^b 


Vhr Vhor 2 (y/ho)^ 


(24) 


Considering H as a function of b one can expand it for small values of the parameter b (till 
linear order) as follows 


H{r,b) ^ Ho{r) + b 


dH 

= Ho{r) + b6H 


6=0 


(25) 


Substituting (24) and (25) into Eq. (21) and keeping only the first order terms in b we arrive 
at the following equation for 6H 


d6H ^ I kr 

oii = —^ _' „ Hq 


dr y/h^r 2 

which can be immediately integrated to obtain 


(26) 


5H = Ha 


1 


kr 


.2 7 {Vhof 

Now we can write the final expresion for H as 


dr 


(27) 


H = Ho{l + bI) 


(28) 


where we denoted by I the following quantity 
^ I f kr ^ Ik 

/ — _ / _ _ _ Wt’ —_ 

2 7 (V^)3 2r-2M 

+ 15M^(r — 


ry/ho{r^ + 5Mr — 30M^) 
2M) ln(2rv^ + 2r - 2M) 


(29) 
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The accuracy of such solutions can be checked by comparison with a numerical one, provided 
that one imposes the same initial condition: Hnumericai{\{rb + ^’c)) = + Vc))- The 

corresponding plots are collected in Fig. 1 for two different values of Ar^. 




Figure 1: In the left panel ~ 4 and Vc ~ 136.52 while in the right panel ~ 10 and 
rc « 542.65. 

Inserting eq.(28) into eq.(20) we obtain a first order nonhomogeneous eqnation for 

+ = + (30) 

dr hr h 

The solution to the homogeneous part of eq.(30) can be obtained analytically by expanding 
around ho for small values of b, in a similar way as it was done above for H. This provides 

= (31) 


We find out a particular solution to the nonhomogeneons equation using the method of 
variable coefficient, i.e. by assuming that F^^^^ = C{r)F^^^, where C{r) is found to be given 
by the following integral 


C{r) 


n + hi\fi \\ r^dh 

\l + V^) \l-bl) 2M-2A/3r^ 


(32) 


Thns the general solntion of Eq.(19) reads now as 


F+ = (A2 + B2C(r))F+„ 


(33) 


where A 2 and B 2 are the two integration constants. In the next section we will match these 
solutions two by two in order to find an analytical expression for the greybody factors that 
characterize the emission of fermions by Schwarzschild-de Sitter black holes. 


4 Greybody factors and energy spectrum 

Let us start by making the observation that the matching of the solutions works only in the 
low-energy regime for which we can expand for small valnes of e the near horizon solutions given 
in (17) and (18) to obtain 

F^ « A^^{1 — ieplnh + ...) (34) 
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F+ « - ieq\nh +...)+ A°^\l + ieqliih + ...) 


(35) 


In order to match with , respectively Ff with F^ we first need to find out the 
shape of the function Ff in the vicinity of the two horizons. The r ^ ri limit of Ff when 
approaches the two horizons is given by the following formula 


lim Ft = Oi {A 2 + Pi B 2 In h) 

r^Ti 


where the constants a* and Pi are defined as 


ai = H^\l-bI) 


r=ri 


Pr = 


1 — y/hp 
1 + 


-2k 


1 + 6 / 


1 - 6/y 2M - 2A/3r3 


r=ri 


Thus, the r ^ rh limit of (33) can be expressed as 


(36) 


(37) 


F^ « Oft A 2 + abPb B 2 In h 

Comparing (34) with (38) will give us the following matching conditions 

A2 = —A^^ B2 = -^A^^ 

Taking now the r ^ limit of (33) yields 

Ff « ofc A 2 + ctcPc B 2 In h 


(38) 


(39) 


(40) 


and comparting this time with the solution at the cosmological horizon given in (35) will lead 
us in the end at 


Taking into account the form of the asymptotic solutions (17)-(18) and the conservation of 
flux we can now compute using relations (39) and (41) the greybody factors for which we found 
the following expression 


r,(e) = l- 


^ f pPc-qPb V 
A™ \pPc + qPb) 


(42) 


The above form of the greybody factors looks very similar with the one obtained in [13, 16] 
for the case of scalar particle emission, for which the authors have also found T^ to be constant 
at low-energies. We see that this feature remains valid also in the case of low-energy fermion 
emission by a SdS black hole. Moreover for fermions the greybody factor (or equivalently the 
absorbtion probability) has a constant value for each mode at very low-energies, for which the 
above (42) formula is valid. However, for fermions Tj has a much more complicated dependence 

8 












on rh and Tc compared to the scalar case for which[13, 16] 


(r-g - rir 
(r^c + rlY 


(43) 


Let us mention that both expressions (42) and (43) for the absorption probability take 
the value 0 in the limit A i—)• 0 (<t4> 6 e-)• 0). This is due to the fact that in this limit ri, i—>• 
2M, Tc e-)• c)o in such a way that brf = ri — 2M. Therefore, a ratio ^ i—)• 0 in (43) and one 
gets expected result. Similarly, the second term in (42) can be expressed as a function of ratios 
First we have to calculate the limits (cf. (29)): limf,^o (^-^c) = ^k, limb_>,o (6/fe) = 0. 

These imply limb_).o i-Ifc ’ (^) ~ ^ which provide the required answer 

limA^o (rj) = 0. 

In Table 4 we present some numerical values of the greybody factors evaluated for the first 
three mods (j = 1/2,3/2, 5/2) using different values for Arg. We observe from the table that for 
each mode the value of Tj becomes higher as we increase the value of the cosmological constant 
A. Another conclusion that can be drawn is the fact that the contribution of the lowest mode 
j = 1/2 to the emission spectra will be the dominant one. 


-wr 

0.001 

0.0025 

0.005 

0.0075 

0.01 


1.14 

2.95 

6.11 

9.42 

12.85 


1.38 • 10"^ 

9.07 • 10-5 

3.84 • 10-"^ 

9.03 • 10-^ 

1.67-10-8 


1.2 • 10"^° 

1.94- 10-9 

1.56 • 10-® 

5.9 • 10-8 

1.2 • 10-^ 


Table 1: The greybody factors for the first three modes (all the numerical values of rj(e) have 
been multiplied by a factor of lO'^). 


Comparing numerically the greybody factors (or equivalently the absorbtion cross section) 
given by eq. (42) and (43) in the massless limit for the lowest angular quantum numbers {j = 1/2 
for fermions, respectively s = 0 for scalars) we obtain that their ratio is approximatively 



T 


s=0 


oc 


—obs 

J 2 
^ahs 
^s =0 


1 

12 


(44) 


For pure Schwarzschild black holes this ratio was showed in ref. [10] to be equal to 1/8. 
Thus we see that in the presence of a cosmological constant a SdS black hole will emit even 
more low-energy scalar quanta as Hawking radiation compared with the radiation emitted by a 
pure Schwarzschild black hole. 

Let us now address the problem of energy spectrum. It is known that the energy emission 
rate for spin 1/2 particles for a SdS black hole has the usual expression[14] 


dE 

de 


1 

27rv‘^ 


+ 1 ) 


exp(e/rff) + 1 


(45) 


where v represents the speed of the fermion and Th is the Bousso-Hawking temperature [38] of 
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the black hole 


1 1 


1 

n 


(46) 


Th 


/i(ro) 47r 



with ro = (3M/A)^/^ obtained from the zero point of the hrst derivative of h{r) . 

The energy spectrnm obtained nsing (45) and (42) is displayed in Fig. 2. This spectra shonld 
be trnsted for qnantitative resnlts only in the low-energy regime. In the intermediate and high- 
energy regimes the spectra is expected to suffer modifications dne to the fact that in these 
regimes the greybody factor will eventually develop a dependence on energy. 

We observe an enhancement of the spectrum with the increasing value of the cosmological 
constant. In the case of asymptotically flat Schwarzschild black holes the energy emission rate 
for fermions vanishes in the limit e —)• 0. We see that this remains valid also for SdS black boles. 
However, for massless scalar particles it turns out that this limit has a non-vanishing value for 
SdS as reported in ref. [16, 17, 18, 14]. Taking the limit of (45) as the energy goes to zero one 
can see that as long as the greybody factor has a constant value the energy emission rate will 
always go to zero regardless of the fermion mass. In the case of massive scalars it was showed 
in ref. [18] that the energy emission rate again goes to zero for a SdS black hole. 



Figure 2: The fermion differential energy emission rate for different values of Ar^. For the left 
panel we have set = 1, respectively = 5 for the right panel. 

We should point out that our analytical results obtained here for low-energies are in good 
agrement with the results reported in [14] where the authors have numerically investigated the 
problem for all types of helds (scalars, fermions, bosons and gravitons). 

5 Conclusions 

In this paper we calculated for the first time, analytically, the low-energy greybody factors for 
fermions emitted by a Schwarzschild-de Sitter black hole with a small positive cosmological 
constant A. We have found that the fermion greybody factors can be approximated by constant 
quantities at low energies. This result is similar to the one derived for scalar particles [13, 16]. 
Our analytical results presented here seems to confirm the numerical studies reported in Ref. [14] 
for very small energies. Another important aspect encapsulated in our results is the fact that 
the ratio fermions to scalars (which is approximatively 1/12) emitted by a SdS black hole is 
smaller than the same ratio[10] (1/8) in the case of a Schwarzschild black hole. 

This work can be continued in a number of ways such as: taking into account the next 
order terms in solving the equation in the intermediate region presented in section 3.2 and 
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thus possible obtaining the energy dependence of the greybody factors; Ending more accurate 
solutions in the regions near the black hole and cosmological horizons and thus refining the 
results presented here. 
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